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Abstract 

Supernova neutrinos, which arrive at Earth earher than hght, allow for the earliest 
determination of the direction of the supernova. The topic of this paper is to study how 
accurately we can determine the supernova direction. We simulate supernova neutrino 
events at the SuperKamiokande detector, using a reahstic supernova model and several 
realistic neutrino oscillation models. With the results of our simulation, we can restrict 
the supernova direction to be within a circle of radius 9°. In several neutrino oscillation 
models, this accuracy is increased to 8°. We also discuss the influence of an accident 
that occurred at the SuperKamiokande detector. After repair of the detector, using 
the remaining PMTs, the accuracy becomes about 12° for no oscillation. 
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1. Introduction 



In order to obtain information concerning core-collapse supernova explosions, which are 
rare in the galaxy, it is very important to observe the light curve of its early phase. It is 
expected that a core-collapse supernova explosion in our galaxy will be detected in the future 
by several neutrino detectors around the world. When a supernova collapse occurs, neutrinos 
produced in the core can escape immediately, because of their very weak interaction with 
matter. Contrastingly, photons do not escape until the shock wave travels from the core 
through the stellar envelope and breaks out of the stellar photosphere. For this reason, 
we can receive the neutrino signal several hours earlier than the light signal. (Of course 
this time delay depends on the size of the envelope. Reference 2) contains a simple model 
of this delay.) In addition, because electromagnetic signals are obscured by dust in the 
interstellar space, it is plausible that some supernova explosions cannot be detected without 
the neutrino signal. Therefore, if we could determine the direction of a supernova explosion 
by its neutrinos, many astronomical observations of its early state would be possible. In 
fact, a world-wide early supernova alert project is running (SNEWS, or SuperNova Early 
Warning System) . 

This problem (of determining the direction to a supernova by its neutrinos) has been 
studied previously in a general manner ^) " ^) . There are two methods to deal with this prob- 
lem. The first method uses the angular distributions of the neutrino reaction products, which 
can be correlated with the supernova direction. Among past works following this approach, 
Beacom and Vogel in Ref. 5) evaluated the numerical integral to find the centroid of the 
Gaussian peak of the reaction products' distribution, with a known fiat background. Their 
result is 69 ~ 5° for SuperKamiokande, and 59 ~ 20° for SNO. The second method is based 
on triangulation using two or more widely-separated detectors. However, this technique was 
shown to be very crude in Ref. 5) {S cos 9 ~ 0.5 for SuperKamiokande and SNO), in contrast 
to previous optimistic estimates given in, for example, Ref. 6). 

In this paper, we simulate a supernova explosion in the galactic plane {D = 10 kpc as 
used in Ref. 5), for comparison) and discuss the statistical error in its direction using the 
first method mentioned above. This approach, which is based on a more concrete simulation 
comparing to the past works, is very precise. In addition, we also consider the more realistic 
case, that of "neutrino oscillation," which is supported by solar " and atmospheric 
neutrino data. This realistic case was not considered in previous works. We expect that 
in the case of neutrino oscillation, the accuracy of the determined supernova direction will 
be better, since the energy of u^s, which give the largest contribution in determining the 
direction, is increased due to the conversion 1/^ i/^,-. (Because i^n^s experience only 
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neutral-current reactions in the supernova and interact with matter more weakly than z/g's, 
z/^^r's reach equilibrium deeper in the core than z/g's, and their temperatures are higher than 
those of z/g's.) 

The construction of this paper is as follows. In §2, a realistic supernova model and 
the neutrino oscillation models used in the simulation are presented. Reactions at the 
SuperKamiokande detector and these cross sections are discussed in §3. Based on these 
models and reactions, we simulated a supernova explosion in our galaxy and events at Su- 
perKamiokande. The results of these simulations are reported in §4. The details of this 
simulation and its results concerning the accuracy of the direction are also reported in this 
section. Finally we discuss our results in §5. 



§2. Supernova model and oscillation parameters 

We use a realistic model of a collapse-driven supernova proposed by the Lawrence Liv- 
ermore group. Time-integrated energy spectra for the case of no oscillation is shown in 
Fig. 1. (See Ref. 15) for details.) In this model, however, the radiation of neutrinos is 
isotropic, despite the fact that supernova progenitors are rotating (which many observations 
have indicated), and they radiate anisotropic neutrinos. For now, we do not consider the 
effect of rotation for simplicity. 

We use four models of neutrino mixing parameters used in Ref. 17) (see Table I), and the 
results of the supernova neutrino oscillation given in Ref. 17) are also used. These models are 
made to agree with the results of the solar and atmospheric neutrinos. " The expressions 
"LMA" and "SMA" indicate the MSW solution of the solar neutrino problem. A recent SNO 
observation, along with other observations, shows that the LMA solution is better than 
the SMA solution, and the SMA solution is not allowed below the 3(7 level (see Ref. 18) and 
references therein). However, we also consider the SMA solution for comparison. The suffices 
"-L" and "-S" attached to "LMA" and "SMA" indicate whether 613 is large or small, where 
^13 is one of three mixing angles of the neutrino mixing matrix. A large (small) ^13 means 
that the "higher resonance" is adiabatic (nonadiabatic) . An adiabatic higher resonance 
enhances the energy of the electron neutrinos, and enhances the event rate of i/e scattering, 
which appears to be strongly correlated with the supernova direction (see below). For a 
review of the MSW effect, see Ref. 20). 
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§3. Expected events at SuperKamiokande 

SuperKamiokande (SK) is a water Cherenkov detector with 32,000 tons of pure water 
located at Kamioka, Japan. The relevant interactions of neutrinos with water are as follows: 

Ue+P^n + e+, (CC) (3-1) 
Ue + e- ^Ue + e', (CC and NC) (3-2) 
Pe + e~^Pe + e~, (CC and NC) (3-3) 

+ e- ^ + e", (NC) (34) 

Ue + O ^ F + e-, (CC) (3-5) 
i?e + ^ N + e+, (CC) (3-6) 

where CC and NC stand for charged current and neutral current interactions, respectively. 

The efficiency of the SK detector is 100% for an electron whose energy is above 5 MeV 
and 50% at 4.2 MeV. (For the energy and angular resolution of SK, we refer to Ref. 21).) 
The energy resolution is ~ 15% for an electron with energy 10 MeV. We display the angular 
resolution as a function of the recoil electron energy in Fig. 2. In this figure, we fitted 
experimental data to the function 83°£'~°-^, where is measured in MeV, which is used in 
our simulation. 

The differential cross section of the electron scattering (3-2)-(3-4) (for its derivation, see 
Ref. 22)), is given by 

da G|r2(l + 2me/Te)3/2 



dcosO 2% 1 + Trie/ El, 



X 



EJ El 



(3-7) 



cos^ = _— — , (3-8) 



E, \T^ + 2me 

where Gp is the Fermi constant, rrie is the electron mass, Tg is the electron kinetic energy. 
El, is the neutrino energy, and the coefficients A, B, and C are given in Table II. 9 is the 
angle between the injected neutrino and the recoil electron (positron). This differential cross 
section is highly peaked in the forward direction, as shown in Fig. 3. With a threshold 
energy (5 MeV), this forward peak is enhanced. 

The differential cross section of the v^P CC reaction (3-1) (see Ref. 23)) is 



da (Jo 



(/^ + 3/) + (/^-/)cose-^ 



d cos 9 2 

r^2{f + h)g[{2Ef^+A){l-cos9)] 



Ef)Ef\ (3-9) 
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+ if + g'mi + cos9)] 

+ if + 3/)[3(Ef + A){1- cos 9) - A] 

+ if - 5')[3(£^i°^ + - cos^) - A] cos 9, 



(3-10) 



where f — l,g — 1.26, and f2 — l^'p — l^n — 3.706 (with fip, /x„ the magnetic moments of the 
proton and a neutron, respectively, in units of the nuclear magneton), M is the nucleon mass, 
and A is the mass difference between a neutron and a proton. The normalizing constant ctq, 
including the energy-independent inner radiative corrections, is 

-o = 5^^(l + 4lJ. (3-11) 

where 9c is the Cabibbo angle (cos^c" = 0.974) and A^^^^ ~ 0.024. This cross section 
(3-9) is expressed to first order in Ey/M, and depends on the zeroth-order positron energy 
Ef^ = Ey — A. As shown in Fig. 4, the differential cross section of the v^P reaction is almost 
isotropic, and the v^P reaction has the largest contribution to the detected events at SK 
[e.g., at Ey = 10 McV, (j{pcp) — 100o"(z/ee )]• For this reason, it is not easy to see the peak 
position of events, information about which enables us to determine the supernova direction 
easily. 

The differential cross sections of the reactions with oxygen (3-5) and (3-6) are unclear, 
because of the uncertainty of the the nuclear part. These reactions are important if we 
consider the case with "neutrino oscillation," because the oscillation enhances the energies 
of and i^e, which contribute to the reactions with the oxygen, and therefore enhances the 
cross sections of these reactions. The differential cross sections with oxygen calculated in Ref. 
24) are highly peaked in the backward direction when the neutrino energy becomes larger, 
and therefore the reactions with oxygen seem to be useful for determining the supernova 
direction. However, the backward peaks of these reactions are not as sharp as the forward 
peaks of the electron scattering reactions. For this reason, we do not consider the reactions 
with oxygen, because we believe that these reactions influence the result only slightly. 

The number of events at SK is calculated in Ref. 17) in the case of a supernova explosion 
at = 10 kpc. We give these results in Table III. (Note that the reactions with oxygen are 
also shown.) 

We cannot expect better accuracy for the direction from SNO observation than from SK 
observation (see Ref. 5), for example), and therefore we do not discuss events at SNO. 
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4. Simulation and results 



We assume that a supernova explosion occurred at D = 10 kpc. To obtain information 
about the direction of the supernova, we have to set the coordinates at SK. Here, we set the 
2;-axis in the upward direction, and we use spherical coordinates {9, (f)) as follows to determine 
the direction: 



We also assume here that the direction of the supernova is {6, 0) = (135°, 270°). In this case, 
the peak position should be (6*, (p) = (45°, 90°). (As we see below, this assumption does not 
influence our flnal results.) 

Recoil electrons and positrons are distributed according to the equations (3-7) and (3-9). 
However, they also experience multiple scattering in water, and we understand this effect as 
follows. An electron or a positron of energy Eg that scatters is characterized by a Gaussian 
probability distribution, whose center value is the original direction of that electron and 
whose one sigma error is 83° x (Ee/MeV)^^^'^ (as shown in Fig. 2). In addition, another 
factor that effects SK events is the energy resolution of the detector. However, we do not 
have to consider the energy resolution, because it represents the accuracy with which the 
detector can determine the electron energies, and in our simulation we do not investigate 
the electron energy distribution. 

First, we consider the case of "no oscillation." We display in Fig. 5 the distribution of 
the events for this case. There seems to be an obscure peak around {6,(f)) = (45°, 90°). To 
analyze this result systematically, wc divide the cos^^ direction into 20 bins and sum up in 
each bin. Figure 6 is the result of this operation. The dashed curve in Fig. 6 corresponds to 
the = 0° peaked case, which can easily be calculated theoretically with cross sections. Then, 
we rotate the coordinates so that the events are seen as in the case that the peak is located 
at = 0° using a least-square method. The best-fltted result is shown in Fig. 7. This result 
is obtained at angles (^^bcstfit, ^bcstfit) = (47.7°, 80.0°). This result is near the "true" values, 
(45°, 90°), but this simulation alone is not sufficient to estimate the errors. We, carried out 
the simulation 1,000 times under the same conditions and obtained a distribution of the 
best-fitted coordinate rotation angles. We denote by ^sn the angle between these best-fitted 
points (6'bestfit; 0bestfit) and the "true" point (45°, 90°), and we measure these discrepancies 
^SN- [We should note that when we analyze using 6'sn, the information regarding the original 
peak position, or (45°, 90°), is lost, and our analysis from this point does not depend on 



X = cos 6 sin 0, 
y — sin 9 sin 0, 



z — cos 9. 



(4-1) 
(4-2) 
(4-3) 
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where the peak position is.] The sohd histogram in Fig. 8 is the distribution of the values 
cos^^SN, and wc fitted this resuh to a Gaussian (dashed curve). The one sigma error of 
this Gaussian is 5cos^^sn = 0.0130. By using the simple formula 5cos^sn — (^^sn/^? we 
obtain dOsti = 9.2°. Then in this method we can determine the supernova direction with an 
accuracy of 9.2°. 

We adapted the same simulations to the other four "neutrino oscillation" cases. We 
summarize these results in Table IV. 

§5. Discussion 

As shown in Table IV, we can determine the supernova direction within ~ 10°. In the 
five models with which we have dealt in this paper, "LMA-L" and "SMA-L" are the best, 
because the electron scattering events are more prominent for them, as shown in Table III. 

Now we give further discussion of three points below, where we assume the case of "no 
oscillation." 

5.1. Oxygen events 

As discussed above and shown in Table III, we must consider the oxygen events, especially 
in the case of "neutrino oscillation." We do not expect that these effects will enable us to 
determine the direction more accurately than in the no oscillation case. Rather, because 
their reaction cross sections are unclear, they may complicate our analyses by acting as 
noise. 

If we would like to determine the supernova direction without the oxygen, whose effect 
is very complicated, one method is to include an "energy cut-off." For example, we only 
consider electrons whose energies are less than 15 MeV, so that oxygen events should have 
little contribution. (However, we cannot be sure that an energy cut-off of 15 MeV is sufficient 
to accomplish this.) With this energy cut-off, we carried out simulations in the same manner. 
For the result obtained in this case, the accuracy is rather worse, S9sn ~ 13°. This is because 
lower energy electrons (positrons) are scattered in the detector, and this makes the angular 
resolution of these electrons (positrons) worse, as shown in Fig. 2. 

5.2. Dependence on distance 

The simulation and its results given above are based on the assumption that the supernova 
exploded at a distance D — 10 kpc. Since the event number at the earth falls off as the 
distance D squared, N oc and the accuracy is proportional to 1/\/N, we expect 
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that the accuracy is proportional to the distance, or, 



5^SN = 9.2° 



( 



10 kpc 



D 



) 



(5-1) 



Actually, wc carried out simulations assuming other distance {D = 5, 7.5, 12.5, 15, 17.5, 20 
kpc) and obtained approximately the same results as above equation. 

5.3. After the accident of SuperKamiokande 

On November 12th, 2001, an accident occurred at the SuperKamiokande detector, and a 
significant part (about 60% of the PMTs) of the detector was damaged. The Superkamiokande 
detector is expected to restart within a year or so, by rearranging the PMTs, whose number 
density is reduced by a half, using the existing resources, and using the same volume of 
water as before. After this repair, the effect of the accident on its performance is expected 
not to be serious for supernova neutrinos, because the fiducial volume will not change, and 
the threshold energy change (from 5 McV to about 7-8 MeV) infiuences the event number 
very little, since the low energy event number is small. (For example, the total event num- 
ber of electron scattering and the inverse f3 decay reaction is 8,318 for a threshold of 5 MeV 
and is 8,165 for a threshold of 8 MeV, based on the calculation of Ref. 17).) The energy 
resolution will be rather large, but this is not a serious problem, since it is proportional to 
square root of the number of PMTs. 

For analysis of the direction, the detector performance is not changed significantly by the 
change in the detector. The reason for this is as follows. First, since the fiducial volume does 
not change, the total event number decreases, only because of the higher energy threshold. 
This effect is very small, as shown above. Second, the angular resolution of the detector, 
which is most important in determining the accuracy, is affected little because the higher 
energy threshold cuts low energy electrons (positrons), which mainly obscure the peak. For 
these reasons, the angular resolution is not strongly affected, and therefore neither is the 
accuracy in determining the supernova direction. 

Actually, our simulation with threshold 8 MeV gives 59sn — 11.9°. This difference comes 
mainly from the relative smallness of the electron scattering events. But this effect cancels 
with the higher angular resolution, discussed as the second reason above, and results in the 
small difference 2.7°. 

We expect that the SuperKamiokande detector will be repaired soon. 
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Fig. 1. Energy spectra of Ug, Ug and of the numerical supernova model used in this paper (, 
where Vx means i^ju^,- and i'ix^t)- The solid, dashed and dotted lines are the spectrum of Vg, 
and Ux, respectively. These spectra are assumed "no oscillation" 
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. 3. Cross section for i/g + e~ — > + e~, with Ei, = 10 MeV. The soHd curve represents the 
cross section without an energy threshold. The dashed curve represents that with an energy 
threshold of 5 MeV, and that threshold is of SuperKamiokande detector. With the threshold, 
only the event of cos 6 > 0.95 can be seen. 
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weak forward peak can be seen, but this peak is not as sharp as scattering events. 
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Fig. 5. Events at the SuperKamiokande detector. An obscure peak around (^,0) ^ (45°, 90°) can 
be seen. 
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Fig. 6. Analyzed figure of the event distribution (Fig. 5). The data points represent the numbers 
of events in each bin, with 20 bins along the cos direction. The dashed histogram corresponds 
to the ^ = 0° peaked case, which can easily be calculated theoretically with cross sections. 
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Fig. 7. The analyzed figure after the coordinate rotation. The rotation angle is determined by 
a least-square method, so that the data points in Fig. 6 can be seen as in the case that the 
peak is located at 9 = 0°. This best-fitted result is obtained at the angles (^bestfit; ^bestfit) = 
(47.7°, 80.0°). 
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Table I. Sets of mixing parameters used in the calculation. Here, dij represents the mixing angle 
of the neutrino mixing matrix and Amfj represents the squared mass difference between the 
i-ih. and j-th. mass eigenstate of the neutrinos. 



model 


sin^ 2^12 


sin^ 2^23 


sin^ 2^13 






Uq problem 


LMA-L 


0.87 


1.0 


0.043 


7.0 X 10--' 


3.2 X 10--^ 


LMA 


LMA-S 


0.87 


1.0 


1.0 X 10-^ 


7.0 X 10-5 


3.2 X 10-3 


LMA 


SMA-L 


5.0 X 10"^ 


1.0 


0.043 


6.0 X 10-6 


3.2 X 10-3 


SMA 


SMA-S 


5.0 X 10"^ 


1.0 


1.0 X 10-6 


6.0 X 10-6 


3.2 X 10-3 


SMA 



Table II. Coefficients for the cross section of ue ue . Here, gv = 2 sin^ 9\y — ^, and gA = 
where 9w is the Weinberg angle (sin^ 9w = 0.23). 
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Table III. Number of events at SuperKamiokande. 



model 


LMA-L 


LMA-S 


SMA-L 


SMA-S 


no osc 




9459 


9427 


8101 


7967 


8036 


Veer 


186 


115 


189 


131 


132 


Z/eC- 


46 


46 


41 


42 


42 




25 


26 


25 


30 


30 




24 


23 


24 


24 


24 




25 


26 


25 


30 


30 


z/^e- 


24 


23 


24 


24 


24 


Oz/e 


297 


214 


297 


108 


31 


Oz/e 


160 


158 


95 


92 


92 


total 


10245 


10114 


8822 


8447 


8441 
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Table IV. Errors of best-fitted angles of the coordinate rotation in degrees. 



model 


66 SN 


LMA-L 


8.1 


LMA-S 


10.7 


SMA-L 


8.1 


SMA-S 


9.3 


no osci 


9.2 
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